Universal Elasticity and Fluctuations of Nematic Gels 



O 

o 

(N 
00 



o 

CO 



T3 

o 
o 



> 

(N 
(N 

o 

(N 
O 
-i— > 



i 

C 

o 
o 



X 
S3 



Xiangjun Xing and Leo Radzihovsky 
Department of Physics, University of Colorado, Boulder, CO 80309 
(Dated: February 1, 2008) 

We study elasticity of spontaneously orientationally-ordered amorphous solids, characterized by 
a vanishing transverse shear modulus, as realized for example by nematic elastomers and gels. 
We show that local heterogeneities and elastic nonlinearities conspire to lead to anomalous nonlocal 
universal elasticity controlled by a nontrivial infrared fixed point. Namely, at long scales, such solids 
are characterized by universal shear and bending moduli that, respectively, vanish and diverge at 
long scales, are universally incompressible and exhibit a universal negative Poisson ratio and a 
non-Hookean elasticity down to arbitrarily low strains. Based on expansion about five dimensions, 
we argue that the nematic order is stable to thermal fluctuation and local heterogeneities down to 
di c < 3. 



Liquid crystal elastomers and gels weakly 
crosslinked networks of liquid crystal polymers — com- 
bine the electro-optic response and thermodynamic phase 
behavior of liquid crystals with the mechanical advan- 
tages of solids, such as rubber and plastics. They there- 
fore hold considerable technological potential. These ma- 
terials exhibit rich interplay between orientational order 
and network elasticity that leads to many unusual prop- 
erties not found in conventional liquid crystals or in con- 
ventional rubber Q. The most striking of these is the 
vanishing stress, cry in response to a finite strain, Uij 
applied transversely to the spontaneous || uniaxial dis- 
tortion that develops below the isotropic-nematic (IN) 
transition. Much progress had been made in understand- 
ing these materials both from the neo-classical theory of 
rubber jjj and from more general symmetry-based elastic 
formulation j^, Q). Many of the intriguing properties of 
nematic elastomers can be traced back to the existence 
of novel nemato-elastic Goldstone mode|5|, |(|, associated 
with a spontaneous breaking of rotational symmetry J3], 
of the amorphous polymer matrix. Although consider- 
able progress had been madc(l), most of the analyses 
had been limited to mean-field treatments. We have re- 
cently developed a fully nonlinear elastic theory of ne- 
matic elastomers 0, that, for one, allows us to assess the 
effects of thermal fluctuations [||. Furthermore, while sta- 
tistically homogeneous and isotropic]^], elastomers are 
locally quite heterogeneous |J. It is essential to study 
the role that network heterogeneity plays in determining 
macroscopic properties of liquid crystal elastomers, and 
this is the goal of the present Letter. 

We find that on scales longer than ^ Z NL ~ if 2 /A, 
(with K an effective Frank nematic modulus and A a 
measure of heterogeneity) even arbitrarily weak hetero- 
geneity qualitatively modifies liquid crystal and elastic 
properties of a nematic elastomer relative to those of the 
ideal homogeneous and isotropic one|| 0|. In particu- 
lar we find that macroscopically a nematic elastomer ex- 
hibits a nonlocal elasticity, characterized by shear mod- 
uli, that vanish as universal power-laws of system size, 
implying a host of exotic elastic behavior: (1) a non- 



Hookean elasticity o~ zz ~ {u zz ) S (with S > 1 universal but 
geometry-dependent), that extends down to arbitrarily 
weak strains, even along the uniaxial axis z[^|, contrast- 
ing strongly and qualitatively with the linear longitudinal 
stress-strain relation of idealized heterogeneity-free elas- 
tomers; (2) a crossover to linear stress-strain relation for 
strain smaller than a critical strain u c zz ~ H 2 / s set by 
an external aligning (e.g., magnetic) H field; (3) a ne- 
matic Frank modulus that diverges as a power-law of the 
minimum of lengths associated with system size, stress 
or aligning field; (4) universal ratios of elastic moduli, 
with a negative Poisson ratio in the plane transverse to 
the nematic axis z, and (5) a macroscopic incompress- 
ibility, independent of microscopic clastic moduli, even 
for elastomers that are compressible on short scales. In 
particular, (4) and (5) predict that a strain u xx , applied 
transversely to the nematic axispO| must be accompa- 
nied by strains u yy = 5u xx /7 and u zz = — \2u xx /7, and 
strain u zz , applied along the uniaxial axis will result in 
Uxx — Uyy = —u zz /2, independent of any other micro- 
scopic details. 

We also find that because of the compliance of the elas- 
tomer network, the nematic order is far more robust to 
the disordering effects of the local random torques by net- 
work's heterogeneities. In fact, a renormalization group 
analysis and expansion about five dimensions suggests 
that orientational order of nematic elastomers is stable 
to weak network heterogeneity |l6| , in strong contrast to 
nematics confined in rigid gels [|l l|, fL2[ 

Although such quenched disorder and thermally-driven 
elastic anomalies, controlled by a nontrivial low tem- 
perature fixed point have been previously predicted in 
a variety of other systems p2[ |lj, pjjj, to our knowl- 
edge, nematic elastomers are a first example of a three- 
dimensional solid, where these exotic theoretical predic- 
tions can be directly experimentally tested. 

As discussed in great detail in Refs. ||, ||, most of the 
properties of nematic elastomers can be captured by a 
purely elastic description in terms of a Lagrangian strain 
tensor u a p = 9 a R-9^R— 5 a p, with the nematic order pa- 
rameter Q a ji = S(n a n/3 — 8 a f}/2>) (S and n, respectively, 
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the magnitude and the unit director characterizing the 
nematic order) integrated out. In this formulation the 
IN transition is signaled by a spontaneous uniaxial dis- 
tortion characterized by a strain oc Q a p, that takes 
place when the effective shear modulus drops below a 
critical value. 

To study the properties of such an oricntationally- 
ordered solid, it is convenient to express its elasticity in 
terms of the strain tensor (and its gradients) measured 
relative to the spontaneously uniaxially strained state. 
The resulting elastic theory resembles that of a conven- 
tional uniaxial solid, but with an essential difference, that 
the shear modulus (A z ±, associated with distortions trans- 
verse to the uniaxial axis identically vanishes || Q. This 
important feature captures, at the harmonic level, the 
soft elastic mode associated with the underlying rota- 
tional invariance of the high temperature isotropic elas- 
tomer solid from which the uniaxial state spontaneously 
emerges. 

Of course real elastomers are only statistically homo- 
geneous and isotropic. As discussed in closely related 
contexts fl2[ |l4| , at long length scales, the elastomer 
network heterogeneities lead to quenched local random 
stress <T Q ^(x) and anisotropy g(x) fields Jl7| 



Ti-R = -u a p(rap(x) ~ (g(x) • n) 



(1) 



that locally distort the elastomer relative to an idealized 
microscopically homogeneous and isotropic state. 

As we will show explicitly below, because of the soft 
mode (vanishing fi z x), such random fields (as well as 
thermal fluctuations) lead to large elastomer distortions 
at which some of the nonlinear elastic terms become com- 
parable to harmonic ones. It is therefore essential to cap- 
ture the full underlying (target and reference spaces [Q) 
rotational invariance (soft mode) of the nematic elas- 
tomer. As we have recently demonstrated, this is quite 
nontrivial and requires a proper treatment of the nonlin- 
ear parts of the Lagrangian strain, as well as keeping of 
the cubic and quartic terms in Lagrange strain tensor^). 
The upshot of that analysis is that for weak heterogene- 
ity, at long length scales elastic and oricntational proper- 
ties of a nematic elastomer |ll| are captured by an elastic 
Hamiltonian 



H = - [B z w 2 zz + Xw% + 2Cw zz wu + 2fiWijWij 



K{V\u z ) 2 } -<t(x)-V_lu, 



(2) 



where 



w zz = d z u z + -(Vim z ) 2 , 

Wij = ^(diUj + djUi) - ^diu z dju z , (3) 

are the rotationally-invariant smectic-like and columnar- 
like nonlinear strain tensors relative to the uniaxial state, 



i and j are summed over indices x and y transverse to 
the uniaxial axis z, and o~i (x) w i(7 i2 (x) + 2gi(x)g z (x), 
that we take to be zero-mean Gaussian quenched random 
field, with spatially independent correlation function 



erj(x)cr,-(x') = A5ij6(x - x') 



(4) 



Our goal here is to study long scale properties of Tt. 
The harmonic correlation functions (w Q (q)u£j(q')) = 
G° /3 (q)(27r) d 5 d (q + q') of the phonon fields u a = (ui 7 u z ) 
can be easily calculated, 



G^(q) - G°g(q) + AqlG° az (q)G%{q) 



(5) 



where G°^(q) are harmonic correlation functions of an 
idealized homogeneous and isotropic elastomer, given by 



G° Z M = 



1 



G° Z M = - 



B z {l-p)ql+Kq±> 
Cq z q, 



(0) 



(A + 2p)q\{B z {\ — p)ql + Kq\) ' (?) 



B z q 2 + Kq\ 



qiqj 



{X + 2p)(B z {l-p)q 2 z +K q \) q\ 



and p = C 2 /B z (\ + 2p) is a dimensionless ratio. In Eq||, 
the second term describes the dominant heterogeneity- 
induced (frozen) deformations, relative to the ideal ne- 
matic elastomer state, with the first giving thermal fluc- 
tuations about this ground state. 

The validity of the perturbation theory in heterogene- 
ity and elastic nonlincarities is controlled by the fluctu- 
ation of the phonon fields. At a harmonic level, at long 
scales, a representative real-space distortion is given by 

i 

KW) = Const.A (^^) 2 L^ d , (9) 

The divergence of these distortions with the size of the 
system Lj_. z for d < 5 signals the breakdown of harmonic 
elasticity on scales longer than nonlinear length scales 
£nl> w hich, in 3d are given by 



(B z (l - /3 ))V4AV2 ' 



¥ NL *K 2 /A. (10) 



To assess the physical consequence of these diver- 
gences, we employ a momentum- shell renormalization 
group (RG), controlled by an expansion in e = 5 — d. 
It is convenient to work with the moduli that renormal- 
ize multiplicatively, 



B = ^(B Z +B ± + 2C), 
ML = \{B Z +B±-2C), 
C = \(B Z -B ± ), 



(11) 
(12) 

(13) 
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where B is the overall bulk modulus, pl the longitudinal 
shear modulus and C the cross coupling. The results 
of the RG coarse-graining procedure are summarized by 
flow equations for the effective elastomer parameters on 
length scale ae l (a a cutoff set by network's mesh size) 

dB 
~dJ 
dC 
~dl 
dfL_ 
dl 
dp 

~dl 

dK 

~dT 
dA 

~dT 

where, 
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(14) 
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Vk 



Vl = ^9l, 



16 5J 



8(1- pl) 9 2 l + 3 (1 + Pi + 2p 2 ^pl) g L g± 



■2 V. 



64 (1 + pi - 2p 2 ^fpl) ql + 96pig± 

gx (4 (2 + 3 pi - 5 pi^fpt) ffL + 5 pi g±) 



4 (2 (l + pi - 2 pi^fpl) g L + 3 pi g x 



The physics is controlled by the flows of two dimension- 
less coupling, gh{p) and g±(l), and two ratios p\(l) and 
pi{l) defined as: 



9l = 



gx 



pi 



P2 



Ap L 

K 3 \ 

gL 
p — j 

PL 
PL_ 

B ' 
C 



K [2B - AC + 3p + 2p L 

i ^ - (15) 

B(3p + 8p L ) + 3p(p L + 2C) - 8C 2 

(16) 
(17) 

(18) 



\ABpI 

Flow equations for pi(l) and p 2 (0 are given by: 
dpi 3 2 

-dT = — 2 ^p^-p^ 

d P2 3 2x 

77" = -l-^^(i-P 2 ), 



(19) 
(20) 



where mechanical stability requires p 2 = C/yJB pl < 1. 
Below 5 dimension, we expect elastic nonlinearities to 
be relevant in the presence of quenched random strains 
and therefore (?l(/) to flow to a positive finite value. This 



together with Eqs|19jj20j implies that at long scales, ratios 
Pl(l) and P2Q) flow to zero, leading to a considerable 
simplification of the flow equations for <?l(Z) (^ 0) and 
9x(l) 



dgL 
dl 

dg± 
dl 



tgL 



zgi 



5.9l (64<7 L 2 + mg L g x + 5lg A 



32 1 



(-64 ffz 



il + %gx) 

752g L g ± 



261g ± 2 ) 



32 {8g L + 3g±) 



(21) 
(22) 



For d = 3 < 5, gz,(0 and g±(l) indeed grow at long 
scales, invalidating harmonic elasticity. Eqs|2l],^, with 
flows displayed in Fig.[l], have four fixed points, Gaus- 
sian (G), Smectic (S), X, and Elastomer (E), that we list 
in Table I. It is clear from Eq.|| that for a vanishing 
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TABLE I: Fixed point couplings and r\ exponents for hetero- 
geneous nematic elastomer. pi i2 = for all fixed points. 

shear modulus p = 0, transverse in-plane phonon modes 
decouple from u z . The remaining longitudinal mode is 
encoded through in-plane density fluctuations wu, that 
can be integrated out, leading to a finite shift in B z , 
with the model reducing to that of a randomly strained 
smectic Q. It is reassuring that this fixed point S, is pre- 
cisely the same as that found in the study of smectics in 
aerogel [O. We are not aware of any physical system that 




FIG. 1: Flow diagram for dimensionless couplings <?l(0 and 



is described by the fixed point X, characterized by gL = 0. 
However, generically elastomers are characterized by a fi- 
nite in-plane shear modulus p and a finite modulus ph 
(gL 7^ 0), and, as can be seen from Fig.[l], S and X are un- 
stable to these couplings. At scales longer than nonlinear 
crossover lengths £atL' the system flows into a globally 
stable zero-temperature fixed point E, controlled by elas- 
tomer heterogeneities and elastic nonlinearities. 

Consequently, at long scales, nematic elastomer can be 
described by an effectively harmonic but nonlocal elastic- 
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ity, with all elastic moduli, except the overall bulk mod- 
ulus B length-scale (wavevector) dependent. Standard 
matching calculation |12] shows that indeed 



K(k) 
M (k) 

C7(k) 
A(k) 
B(k) 



NL\ j 



Ao(fc±^ L )- r,A /A[(^^ L 
Bo, 



NL\ 



(23) 
(24) 
(25) 
(26) 
(27) 
(28) 



where the anisotropy exponent £ = 2 — (77^ + t]k)/2. 
Above scaling functions f a [v] {a = K, fi, /Ul,A) are u- 
independent for large w and scale as v ±ria ^ for a vanish- 
ing w, such that for k^ — ► limit, the elastic moduli are 
fc^-independent. Since this predicts that at long length 
scale, all elastic moduli are vanishingly small compared 
to the wavevector-independent bulk modulus B, nematic 
elastomer is strictly macroscopically incompressible, as 
advertised in the introduction. 

Furthermore, from definitions of the coupling g^ and 
gx we find that the ratio between shear moduli \i and 
\Xh approaches a universal value g]_/ 1 g* L = 24 at the fixed 
point E. This, together with a straightforward analysis 
and Eq.flra predicts that a uniform strain u xx JlOj] leads to 



AH 



' U!) 



2fi R 



12 



(29) 
(30) 



with negative and positive universal Poisson ratios. 

Another fascinating implication of Eqs.|23|-p8| is a 
strictly nonlinear stress-strain response down to an ar- 
bitrarily weak stress a zz < ctnl = ±@- To 
show this we note that a zz cuts off the singular k depen- 
dence (Eqs.|23-27) of elastic moduli on scales longer than 
£ = TO'J'ff-', v = 1/(2 - m ), thereby replacing 
it by a singular er zz dependence ~ (o- zz /o- NL y ia " .jl5| In 

particular we find a non-Hookean response u Z2 ~ o~ X J z , 
1/(5 = 1 — T)l/(2— rjjc), that within (inaccurate) e = 2- 
expansion gives (5 = 157/151 unimpressively close to 1. 

Our results rely on the assumption of a stable long- 
range nematic order fi~6)| , requiring convergent orienta- 
tional fluctuations 



(|H 2 ) * d^l 2 } ~ L^+^ 



(31) 



Thus a necessary condition for stability of nematic order 
is t\k + V± ~ 2 < 0, equivalent to the requirement on 
the anisotropy exponent C > 1. This is satisfied for d > 
di c = 17/56, suggesting a stability of nematic order for 
3d elastomers. 

Experimentally, elastomers crosslinked in the isotropic 
exhibit a polydomain nematic order with a typical micron 



size correlation length||, [T(|. This can be reconciled with 
the above prediction of the e-expansion by heterogeneity 
that is not weak, or by appealing to glassy dynamics that 
prevents a full equilibration on experimentally relevant 
time scales. This hypothesis can tested by cooling in the 
presence of a strong aligning field. 

To sum up, we have studied nematic elastomers, taking 
into account their heterogeneity, and have shown that it 
leads to many striking elastic properties, that should be 
readily experimentally testable. 
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